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Self-Similar Solutions for ADAF with Toroidal Magnetic Fields 
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Abstract 

We examined the effect of toroidal magnetic fields on a viscous gaseous disk around a central object 
under an advection dominated stage. Wc found self-similar solutions for radial infall velocity, rotation 
velocity, sound speed, with additional parameter /? [= c^/(2Cs)], where ca is the Alfvcn speed and Cg is the 
isothermal sound speed. Compared with the non-magnetic case, in general the disk becomes thick due to 
the magnetic pressure, and the radial infall velocity and rotation velocity become fast. In a particular case, 
where the magnetic field is dominant, on the other hand, the disk becomes to be magnetically supported, 
and the nature of the disk is significantly different from that of the weakly magnetized case. 
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1. Introduction 

Accretion-disk models have been extensively studied 
during these three decades (see Kato et al. 1998 for a 
review), and several types of models were proposed. 

One is the classical standard disk, where the accretion 
rate is subcritical, the disk is geometrically thin and op- 
tically thick, Keplerian rotating with negligible infall ve- 
locity, and truncated at the marginally stable orbit. Such 
a standard disk was studied for the Newtonian case (e.g., 
Shakura, Sunyaev 1973) and for the relativistic case (e.g., 
Novikov, Thorne 1973). 

When the accretion rate is so low, the disk is sup- 
posed to be in the optically-thin advection-dominated 
state, where the disk is optically thin with insufficient 
cooling and extends down to the surface of the cen- 
tral object (e.g., Ichimaru 1977; Narayan, Yi 1994). 
Such an optically-thin advection-dominated accretion flow 
(ADAF) has been examined for the Newtonian case (e.g., 
Narayan, Yi 1994, 1995), under the pseudo-Newtonian po- 
tential (e.g., Narayan et al. 1997), and for the Kerr metric 
(Gammie, Popham 1998; Popham, Gammic 1998). 

When the accretion rate is so high, on the other hand, 
the disk is supposed to be in the optically-thick advection- 
dominated state, where the disk is optically thick with 
photon trapping and extends down to the surface of the 
central object (e.g., Abramowicz ct al. 1988; Eggum et 
al. 1988). Such an optically-thick ADAF (slim disk) has 
been examined for the Newtonian case (e.g., Abramowicz 
et al. 1988; Watarai, Fukue 1999; Fukue 2000), under 
the pseudo-Newtonian potential (e.g., Szuszkiewicz et al. 
1996; Artemova et al. 2001), and for the Kerr metric 
(e.g., Abramowicz et al. 1996; Igumenshchev et al. 1998; 
Bcloborodov 1998). 

Disk models with outflows were also proposed. For 
an example, Blandford and Begclman (1999) proposed 
an adiabatic inflow-outflow solution (ADIOS) for an 
optically-thin state, where the wind is driven by the hot 



gas pressure (see also Blandford, Begclman 2004). While, 
for an optically thick state, Fukue (2004) proposed a crit- 
ical accretion disk, where the wind is driven by the radi- 
ation pressure (see also Lipunova 1999; Kitabatake et al. 
2002). 

Also investigated was a convection-dominated accretion 
flow (CDAF), where the disk is convective in the radial 
direction (Narayan et al. 2000; Quataert, Gruzinov 2000). 
In the Newtonian case, self-similar solutions were found 
when convection moves angular momentum inwards. 

The effect of magnetic fields has also been studied, in 
particular, in relation to magneto-rotational instability 
(MRI) (Balbus, Hawley, 1998), and an accretion flow in 
which the magnetic forces dominate over the thermal and 
radiation forces is called a magnetically-dominated accre- 
tion flow (MDAF) (e.g., Meier 2005; see also Shadmchri, 
Khajenabi 2005). However, the intermediate case, where 
the magnetic force is comparable to other forces, has not 
been studied well. Moreover, the mass loss and escape of 
magnetic fields for MDAF were not considered. 

When the sufhcient amount of chaotic magnetic fields 
are created by, e.g., MRI and the ohmic heating operates 
effectively, electrons are heated as well as ions, and ADAF 
conditions would not be realized (e.g., Bisnovatyi-Kogan, 
Lovelace 1997, 2001). We do not consider, however, such 
chaotic fields as well as the ohmic heating in the present 
analysis. Incidentally in this paper, we consider the effect 
of global toroidal fields on the advection-dominated disk, 
taking into accout the mass loss as well as the field escape. 

In the next section the basic equations for the present 
purpose are presented. In sections 3 and 4, self-similar 
solutions are shown. We discuss several properties of the 
present model in section 5. The final section is devoted to 
concluding remarks. 
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2. Basic Equations 

Let us suppose a gaseous disk rotating around a 
Sciiwarzscliild black hole of mass M . The disk is assumed 
to be in an advection-dominated state, where the viscous 
heating is balanced with the advection cooling. In a cylin- 
drical coordinates (r, (/3,z), we vertically integrate the flow 
equations. Furthermore, the flow is assumed to be steady 
and axisymmetric {d/dt = djdip = 0), and all flow vari- 
ables are a function only of r. We ignore the general rela- 
tivistic effects and use Newtonian gravity. We adopt an a 
prescription (e.g., Shakura, Sunyaev 1973). As magnetic 
fields, we consider only toroidal fields B^p. 

For such a disk, the continuity equation with mass loss 

is 

i-^(rEi;,) = 2pi7, (1) 
r dr 

where is the radial infall velocity, p the mass-loss rate 
per unit volume, H the disk half-thickness, and E the sur- 
face density, which is defined as S = 2yoiJ, p being density. 
The radial momentum equation is 
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GM 1 d 



1 d 



^' dr f ^ ^i^^^') f ~2Ed;;^^^^)'^^^ 
where is the rotation velocity, Cg the isothermal sound 
speed, which is defined as Cg =pgas/p, Pgas being the gas 
pressure, and ca is Alfven speed, which is defined 
B'^/{Attp) = 2pmag/p, Pmag being the magnetic pressure. 
On the right-hand side of equation (2), the third term 
implies the pressure gradient force, while the magnetic 
force is designated as the fourth and fifth terms. 

As an alpha prescription, we consider two extreme 
cases: the riy9-component of the viscous stress tensor is 
proportional to the gas pressure or to the total pressure. 
In addition, instead of t^p = rjrdQ/dr = —ap, where 77 is 
the viscosity and a is the viscous parameter (Shakura, 
Sunyaev 1973), we adopt the form of = ^l^^a{p/p), 
where ly is the kinematic viscosity (Narayan, Yi 1994). 
That is, two cases we consider are 
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PmagJ 



case 1, 
case 2. 



(3) 



Hence, the angular momentum transfer equation is 
'aEc^r^ dn" 



riJvr— [rv^) = — 
dr dr 



ilu dr 



case 1, 



(4) 



where (= v^/r) is the angular speed, and JIk (= 
yjGM/r^) the Keple rian angular speed. 

In both cases, the hydrostatic balance in the vertical 
dirction is integrated to 
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Here, we introduce the parameter (3 by 
2 

r),„„„ \ I r \ \ 

/3 = 



Pmag 
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(5) 



(6) 



which is the degree of magnetic pressure to gas pressure. 
We assume this ratio is constant through the disk. 
The energy equation becomes 



7 — 1 dr 

where 7 is the ratio of specific heats. The advection 
parameter / measures the degree to which the flow is 
advection-dominated (Narayan, Yi 1994), and is assumed 
to be constant. 

Finally, since we consider only the toroidal field, the 
induction equation with field escape can be written as 



— (w^B^) =B^, 



(8) 



where B^p is the field escaping/creating rate due to mag- 
netic instability or dynamo effect. This induction equa- 
tion is rewritten as 



dr 



dVr 

dr 



2 ^PH 



(9) 



Self-Similar Solutions for Case 1 



We first consider an advection-dominated accretion 
flow, including toroidal magnetic fields, for the case 1, 
where the viscosity is assumed to be proportional to the 
gas pressure. Under the self-similar treatment by Narayan 
and Yi (1994), the velocities are assumed to be expressed 
as follows, 
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(10) 

(11) 
(12) 

(13) 



where coefficients Ci, C2, and C3 are determined later. 

In addition, the surface density S is assumed to be a 
form of 



E(r) = S]or^ 



(14) 



where Eg and s are constants. Then, in order for the self- 
similar treatment to be valid, the mass-loss rate per unit 
volume and the field escaping rate must have the following 
form, 



(15) 
(16) 



where po and Bq are constants. 

Using these solutions, from the continuity, momentum, 
angular momentum, hydrostatic, energy, and induction 
equations [(1), (2), (4), (7), and (8)], we can obtain the 
following relations. 
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Ci = 3(s+ 1)C3, 



H/t = v/(l + /3)c3, 



3-7 2 _ „ 

7777^1 = 6 ci, 

7-19/ 
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ciaGM, /47ri;o 



/3c3 



(17) 

(18) 
(19) 
(20) 

(21) 
(22) 



As is easily seen from equation (20), in general the disk 
thickness becomes large due to the magnetic pressure for 
the weakly to moderately magnetized cases of /3 ^ 1. In 
addition, for s = —1/2, there is no mass loss, while there 
exists mass loss (wind) for s > —1/2. On the other hand, 
the ecapc and creation of magnetic fields are balanced 
each other for s = 2. 

Considering these properties, we found two types of so- 
lutions under the present treatment. First is the conical 
disk without accretion in the limit of a small alpha, where 
the conical disk is supported by rotation, the gas pressure, 
and the magnetic pressure. In this case, po = and Bq = 0. 

Another type is the advection disk under dynamo action 
in the case of a finite alpha, where the toroical magnetic 
field is injected (and escaping), due to the action of the 



dynamo effect and instability. In this case, p cc r 



and oc r^'^ 



-5)/2 



1-5/2 



in particular, in the case of no wind 

-11/4 



(s = —1/2), p = and B^ cx r' 

Now, the relations (18), (19), and (21) determine the 
constants q. 

When a = at the no- accretion limit, these coefficients 
become 

^ (23) 
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/3 + 3e" 



Be" 



-3e" 
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where 



(l + s)( 



1 - s 



-l3 + 3e" 



23-7 1 



(24) 



(25) 



(26) 



9 7-1./ 

On the other hand, when a ^ 0, these coefficients be 



come 



Cl 




/3 + 3e" +18a2 



-/3 + 3e' 



(27) 



1 



(28) 



C3 



9(l + s)a2 




/3-h3e" -Hl8a2 



1 



/3 + 3e" 



(29) 



The parameters of the model are the ratio of the specific 
heats 7, the standard viscous parameter a, the magnetic 
pressure fraction /?, the advection parameter /, and the 
mass-loss parameter s. In the case without mass loss, 
s = —1/2, as already stated. In addition, the effect of 
mass loss is similar to the advection effect (Fukue 2004). 

Examples of coefficient c^'s are shown in figure 1 as a 
function of the advection factor / for several values of (3. 
Other parameters are 7 = 4/3, a = 1, and s = —1/2 (no 
wind) . 

Figure la shows a coefficient cia, which is a ratio of the 
radial infall velocity to the Keplerian one. Although the 
radial infall velocity is generally slower than the Keplerian 
speed (= GM /r), it becomes large with a or /. In addi- 
tion, the radial infall velocity becomes fast as the toroidal 
magnetic field becomes large. This is because, in the case 
of s = —1/2, the magnetic tension term dominates the 
magnetic pressure term in the radial momentum equation 
(2), that assists the radial infall motion. 

Figure lb shows a coefficient C2, which is a ratio of the 
rotation velocity to the Keplerian one. As is seen in fig- 
ure lb, the rotation velocity decreases as the advection is 
superior in the weakly to moderately magnetized cases. 
This is because the disk gas must rotate faster than the 
case without the magnetic field due to the action of the 
magnetic tension. On the other hand, in the strongly mag- 
netized case (/3 ~ 10), the rotation speed becomes faster 
than the Keplerian speed (super-Keplerian), in order for 
the centrifugal force to be balanced with the magnetic 
force. 

Figure Ic shows a coefficient ^/cz^ which is a ratio of the 
sound speed to the Keplerian one. As the advection pa- 
rameter becomes large, the sound speed as well as the disk 
thickness become large. As for the effect of the toroidal 
magnetic field, the sound speed also increases as the mag- 
netic field becomes large. In addition, as already stated, 
the disk thickness becomes large due to the effect of the 
magnetic pressure. 

4. Self-Similar Solutions for Case 2 

We next examine an advection-dominated accretion 
flow, including toroidal magnetic fields, for the case 2, 
where the viscosity is assumed to be proportional to the 
total pressure. As is easily seen from equation (3), we ob- 
tain the self-similar solutions for case 2, if we replace a in 
equations (4) and (7) by a{l + f3). 

Hence, the solutions at the no-accrction limit of a = 
are 
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C3 



(i + 5)(tt|-^ + 3^" 



[Vol. , 
(30) 

(31) 
(32) 



On the other hand, the sohitions for a 7^ become 
1 



Cl 



3a2(l+/3) 



1-s 

'1-s 



1 + s 



/3 + 3e"j +18a2(l + /3)2 
-/3 + 3e" 



(33) 



3a2(i+/j)2 



1-s 



1 + s 



-/3 + 3e" +18a2(l + /3)2 



1-s 



/3 + 3e" 



(34) 



C3 = 



9(l + s)a2(i + /3)2 



1-s 



1 + s 



-/3 + 3e" +18a2(l + /3)2 



1-s 



-/3 + 3e" 



(35) 



Examples of coefficient c^'s are shown in figure 2 as a 
function of the advection factor / for several values of /3. 
Other parameters are 7 = 4/3, a = 1, and s = —1/2 (no 
wind) . 

In the weakly to moderately magnetized cases, the be- 
havior of the solutions for case 2 is qualitatively similar to 
those for case 1. In the strongly magnetized case, however, 
the behavior is drastically different between two cases, as 
will be discussed in the next section. It should be stressed 
here that both ci and .Jci become small as /3 increased. 

It should be noted that in general case of viscosity: 

^ = ^^K'"Pga«(Pgas+Pmag)'"'', (36) 

where [i is constant, we easily obtain the solutions, if we 
replace a in equations (4) and (7) by a(l + . 

5. Discussion 



Fig. 1. Numerical coefficient c^'s as a function of the advec- 
tion parameter / for several values of /3 in case 1; (a) cia, (b) 
C2, and (c) ,Jc^. The disk density proffie is set to be s = —1/2 
(no wind), the ratio of specific heats is set to bo 7 = 4/3, and 
the viscous parameter is o = 1. The value of /3 are 0, 1, and 
10 from bottom to top of each curve in each panel. 



We briefly note several other properties of the present 
self-similar solutions for ADAF with toroidal magnetic 
fields. 

At first, we show the temperature of the present self- 
similar advection-dominated disk with toroidal magnetic 
fields. 
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Fig. 2. Numerical coefficient c^'s as a function of the advec- 
tion parameter / for several values of /3 in case 2: (a) cia, (b) 
C2, and (c) -y/ca". The disk density proffie is set to be s = —1/2 
(no wind), the ratio of specific heats is set to be 7 = 4/3, and 
the viscous parameter is a = 1. In (a), the value of f3 are 0, 1, 
and 10 from bottom to top of each curve in the panel. In (b) 
and (c), the value of /3 are 0, 1, and 10 from top to bottom of 
each curve in each panel. 



In the optically thin case, where the gas pressure is 
dominant, the isothermal sound speed Cg is expressed as 

GM 



(37) 



where T is the gas temperature, TZ the gas constant, and /2 
the mean molecular weight. In this case, the temperature 
is expressed as 



T- 



- C3; 



cV2 r. 



2.706 X 10 



12 C3 



K, 



(38) 



TZ/^i r ■ r/rg 

where Vg {= 2GM/c^) is the Schwarzschild radius of the 
central object and p, — 0.5. This has a similar form with 
that of the non-magnetic case, but the coefficient C3 im- 
plicitly depends on the magnetic field. That is, the gas 
temperature increases as the magnetic field increases. 

In the optically thick case, where the radiation pressure 
dominant, on the other hand, the sound speed is related 
to the radiation pressure at the disk equator by 

{aT^/3)2H _ ,_GM 

~ "^s ~ „ ; 



where a is the radiation constant, and Tc the central tem- 
perature. Then, in terms of the optical depth r (= KpH), 
the effective temperature Tcff at the disk surface is ex- 
pressed as 



9 1 4 3 
T 4 



C3 



(40) 



where Le (= ^ircGM / k) is the Eddington luminosity of 
the central object. Hence, the effective temperature be- 
comes 



Tcff = 2.181 X 10' 



C3 



M 



1 

-1/4 



1/8 



-1/2 



K. 



(41) 



In this case, the effective temperature very weakly de- 
pends on the strength of the magnetic field, and therefore, 
it may be slightly changed by the existence of toroidal 
magnetic fields. 

The strength of toroidal magnetic fields becomes 



/3C3 



v/(l + /3)c3 



s/2-1 



(42) 



In the case without mass loss (s = —1/2), the magnetic 
field depends on radius as cx r~^^^. This dependence is 
consistent with the usual MDAF (Meier 2005; Shadmehri, 
Khajenabi 2005). 

We briefly discuss the extremely magnetized case of 
large (3 limit. In case 1, where the viscosity is propor- 
tional to the gas pressure, at the limit of large /3, the 
coefficients are approximately expressed as 



Cl 



3a2^' 



(43) 
(44) 
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(45) 



In this case, ci could highly exceed unity, which means 
that the gas dynamically infalls with super-Keplerian 
speed due to the magnetic force. Similary, the gas rotates 
with super-Keplerian speed in this case. 

In case 2, where the viscosity is proportional to the 
total pressure, at the hmit of large f3, the coefhcients are 
approximately expressed as 



Cl 



C3 ' 



Vl + 18a2 + l 



3a2 



£"ci 
ci 



(46) 
(47) 

Hence, in this case, Ci becomes ceiling, while C2 and C3 be- 
come sufficiently small. These properties are qualitatively 
consistent with the results by Oda et al. (2005). 

Finally, in general case, where u cx p^as(Pgas +Pmag)^"^, 
at the limit of large /3 and for < /i < 1, the coefficients 
are approximately expressed as 



Cl ' 



c^ 



C3 ' 



3^ 



1 



(49) 
(50) 
(51) 



9(l + s)a2' 

This means that in the magnetically dominated disk of 
large /3 the nature of the disk is drastically changed at 
H=l/2. For fi>l/2, the gas rotates with super-Keplerian 
speed, and the disk thickness becomes large, while the 
quantities are very small for fi < 1/2. 

Here, we briefly mention on the magneto-rotational in- 
stability (MRI) in relation to the above magnetically sup- 
ported disk. The structure with a strong pure toroidal 
field would be unstable because of a variant of MRI, espe- 
cially, in the interior of stars and accretion disks (Tayler 
1973; Spruit 2002; Akiyama et al. 2003; Ardeljan et al. 
2005; Balbus, Hawley 1992). In the case of accretion disks, 
the growth time of MRI is about where J7 is the an- 
gular speed of the disk, while the inflow time is about 
l/(af2) in the advection dominated accretion flow consid- 
ered in the present paper. Hence, the toroidal magnetic 
field would infall as it grows via MRI, and the strongly 
magnetized disk could be marginally stable. In addition, 
from weakly to moderately magnetized cases, self-similar 
solutions with toroidal fields may be possible because the 
fields linearly grow. 

Furthermore, in numerical simulations of core-collapse 
supernovae with magnetic fields (Akiyama et al. 2003; 
Ardeljan et al. 2005), the linear growth of toroidal fields 
is terminated by the development of MRI, leading to dras- 
tic acceleration in the growth of magnetic energy. In the 
case of magnetized accretion disks, however, a recent re- 
sistive MHD simulation of black hole accretion flows by 
Machida et al. (2005) shows that MRI may be suppressed 



in the strongly magnetized case due to the magnetic ten- 
sion and pressure, and the magnetically supported disk 
may be realized (see also Oda et al. 2005). Hence, we 
consider whether the disk with strong toroidal fields is 
unstable or not is still controversial. 

6. Concluding Remarks 

In this paper we have examined the effect of toroidal 
magnetic fields on advection-dominated accretion flows. 
We found self-similar solutions, using the similarity tech- 
nique in analogy to the self-similar solutions by Narayan 
and Yi (1994). Self-similar solutions would be valid in the 
region from about 10 Schwarzschild radii to about 1000 
Schwarzschild radii, in the case that the total disk size is 
about 10^ Schwarzschild radii. In the present intermedi- 
ate case on the strength of the magnetic field, the disk 
structure, e.g., the coefficient q's, slightly depends on the 
strength of the magnetic field. However, the tendency of 
the effect of the todoidal field would become clear in the 
present study; it raises the infall velocity, the rotation ve- 
locity, and the sound speed in the case that the magnetic 
tension is dominant than the magnetic pressure, and vice 
versa. 

The properties of the present model, e.g., the magnetic 
field, are similar to the results of other MDAF (e.g., Meier 
2005; Shadmehri, Khajenabi 2005; Oda et al. 2005), al- 
though in the present model the toroidal magnetic field is 
dominant. 

This model was thought about a toroidal magnetic field 
in defiance of a relativity effect. If the central object is rel- 
ativistic, the gravitational field should be changed. This 
problem is remained as a future work. 
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